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Abstract 

New scalar structure functions with different sign-symmetry properties are defined. These struc- 
ture functions possess different scaling exponents even when their order is the same. Their scaling 
properties are investigated for second and third orders, using data from high-Reynolds-number 
atmospheric boundary layer. It is only when structure functions with disparate sign-symmetry 
properties are compared can the extended self-similarity detect two different scaling ranges that 
may exist, as in the example of convective turbulence. 

PACS numbers: 47.27.Gs,47.27.Jv,47.27.Nz 
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I. INTRODUCTION 



A problem of broad interest is the advection and diffusion of passive scalars in a turbulent 
flow. The classical paradigm of a passive scalar is the temperature field 0(x, t) when the 
heating is small. The temperature increments A r = #(x + r) — 9{r) have been studied in the 
literature Q, 0, in search of scaling in the intermediate range of scales that are unaffected 
directly by either the stirring mechanism or diffusive and viscous effects. For examples of 
early and recent experimental studies, see |4j and p, respectively. The following two types 
of moments of A r , the so-called structure functions of order n, have been employed: 

S n (r) = ((6>(x + r)-0(x)) B ) = (1) 
S w (r) = (|0(x + r)-0(x)r) = <|A r r>- (2) 

Here, r is the separation vector between two spatial positions, and (■) defines a suitable 
ensemble average. For convenience we will call the normal structure functions and (0) 
the absolute structure functions. When r = |r| is small compared to the large scale L, both 
structure functions are homogeneous, i.e., independent of x. Clearly, (JTJ) and (j2J) coincide for 
even n. However, remembering that normal odd moments are zero in the absence of a mean 
temperature gradient while that is not so for absolute odd moments, one may expect, for 
odd n, that there might be perceptible differences in the two classes of structure functions. 
It is useful to quantify these differences and stress the reasons why they might be important. 

The following two comments put the present work in perspective. The first concerns 
the extraction of scaling exponents of structure functions using the Extended Self-Similarity 
(ESS) Instead of examining the scaling of normal structure functions, S n (r), with respect 
to the scale separation r directly, the practice is to examine the scaling relative to another 
structure function, say S m (r), m ^ n. This usually leads to the extension of a possible 
algebraic scaling range, and the relative scaling exponent, S n (r) ~ S m (rY n ' m where ( njm = 
(n/(m, can be obtained with greater confidence. In the literature, the implementation of 
the method has often mixed up normal structure functions and absolute structure functions 
without exploring the differences between them. Further, it was recently shown that in the 
presence of convection, ESS fails to show the existence of two distinct scaling ranges — the 



0. 



nearly passive behavior at small scales and the dominance of buoyancy at large scales 
This point is illustrated in Fig. In the top figure, we show the second order structure 
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function with clearly separated scaling ranges; in the bottom figure, the corresponding ESS 
plot is shown to result in a line of nearly constant slope, without making the necessary 
distinction between the regions marked A and B. (The second set of data in the bottom 
figure corresponds to another set of measurements, with qualitatively similar conclusion.) 

The second comment is that the structure functions, Eqs. JI} and (J2J), have different 
sign-symmetry properties. One has to distinguish the sign-symmetry with respect to the 
reversal of increments A r — > — A r from that with respect to spatial reflection (also known 
as parity) — that is, for the transformation r — > — r. The absolute structure functions (0) re- 
main unchanged for all orders n under both transformations; they have even sign-symmetry. 
The same is true for even-order normal structure functions but the odd-order normal struc- 
ture functions change sign under both transformations; they have odd sign-symmetry or 
sign-antisymmetry. Note that A r — > — A r does not follow from r — > — r although, for homo- 
geneous turbulence, both transformations have the same effect on (0). We wish to comment 
on ESS in the light of the sign-symmetry properties of structure functions. For this purpose, 
it is convenient to introduce new types of structure functions which explicitly emphasize the 
sign-symmetry with respect to the increment A r . This is our basic goal. 

These new quantities, along with ((TJ and (J2J), will be calculated from temperature data 
from a high-Reynolds-number atmospheric boundary layer. The Taylor microscale Reynolds 
number is about 3500. Since the experiments have been described in some detail in [8], we 
shall mention only a few details here. Measurements were made in a boundary layer above 
salt flats of the Dugway Proving Ground in Utah, at a height of 1.75 m above the ground. 
The ground was smooth on the order of a millimeter. Taylor hypothesis was used and 
Urms/U was about 7% with the mean speed U = 0.72 ms" 1 . Measurements were made at 
various times of the day, covering intensely convective motion in late afternoon, essentially 
neutral conditions in the evening hours and somewhat stable conditions until about 11 PM. 
The data records chosen for analysis corresponded to constant wind conditions in magnitude 
as well as direction. Temperature fluctuations were measured by two cold wires of 0.6 fim 
diameter and 1 mm length. The data acquisition system was a standard constant-current 
anemometer system operated at a small enough current to minimize velocity sensitivity. 

The paper is organized as follows. In Sec. |Hj we will construct the new types of structure 
functions, in line with the approach of for velocity increments. Since these structure 
functions are of intrinsic interest, they are presented first before considering ESS. Section 



3 



IIHI contains the results of the ESS data analysis for structure functions with even and odd 
sign-symmetry. In Sec. IIVI we discuss a possible explanation of the effects observed in terms 
of the SO (3) rotation group decomposition of structure functions of different sign-symmetry. 
Some concluding remarks are presented in Sec. 



II. SIGN-SYMMETRY WITH RESPECT TO THE TEMPERATURE INCRE- 
MENT 

Let the probability density function (PDF) of the temperature increment at fixed distance 
vector r be given by /(A r ). First, we can define ((TJ) and (J2J) as 



S n = / A n r f(A r )dA r (3) 

J — oo 

/oo 
|A r | n /(A r )dA r . (4) 

The PDF can be decomposed into its symmetric and antisymmetric parts with respect to 
the increment as 

WAr) = /(A,) + (5) 

Note that f a does not have the positive-definite property of a PDF. 
We can also define the positive and negative parts of the PDF as 

p(A r ,A r >0) = /(A r ) (7) 
n(A r ,A r >0) = /(-A r ), (8) 

and define moments of A r with respect to f s , f a , p and n, respectively, by the following 
relations: 

/•CO 

P n = A? p(A r )dA r , (9) 
Jo 

/•CO 

N n = A? n(A r )dA r , (10) 
Jo 

/•CO 

S n , s = 2 / A™ f s (A r )dA r , (11) 
Jo 

/•oo 

S n , a = 2 / A"/ a (A r )dA r . (12) 
Jo 
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The following relations are now valid: 



S n>s = P n + N n = S H (13) 

S n ,a = Pn~N n . (14) 

For odd values of n = 2k + 1, Eq. ()14j) reduces to the normal odd-order structure function 

S2k+l,a — P2k+1 — N 2 k+1 = S2k+1 , (15) 

whereas, for even n = 2k, Eq. (J14)) is a new structure function which is even order, but 
sign-antisymmetric. It is not possible to have this combination in terms of either normal 
or absolute structure functions. 

On the basis of the discussion of the sign-symmetries of f s and f a , using the definitions 
(J7J) and (JHJ), it is clear that 

POO 

P n + N n = / A?(p(A r )+n(A r ))dA r 
Jo 

POO 

= / A?(/(A P )+/(-A P ))dA P (16) 
Jo 

is sign- symmetric, while 

POO 

P n ~N n = / A n T (p(A r )-n(A r ))dA r 
Jo 

POO 

= / A?(/(A r ) - /(-A r ))dA r (17) 
Jo 

is sign- antisymmetric. 

To illustrate the differences between normal and absolute structure functions, we plot 
in Fig. |21 the logarithmic local slopes for four types of moments: P n , N n , S njS and S nj0i . 



Figure 2(a) show these quantities for n = 2. Here S^a is the normal (and absolute) second 
order structure function; |S2, a | = N2 — P2 is the newly defined sign- antisymmetric second 
order structure function. Figure p(b)| shows the same four quantities for n = 3. Although 
the scaling is not impeccable even at this high Reynolds number (see Ref. for comments 
in this regard on the scaling of velocity structure functions), a scaling tendency can be 
discerned in the range between a few mm and a few cm. 

Independent of this detail, it is clear that the exponents, if one were to assign nominal 
values in the scaling range, are not the same for all the different structure functions of the 
same order. In particular, the second-order sign-antisymmetric structure function IS2J = 



A^2 — Pi has a substantially larger exponent than the classical exponent of 2/3. While 
5*3,(1 = 5*3 has a scaling exponent close to the Kolmogorov prediction of 1, the other curves 
have measurably smaller scaling exponents. In particular, the absolute structure functions 
have smaller scaling exponents than the normal structure function. 

At the least, these features are disconcerting to anyone interested in scaling exponents. 
An understanding of differences in the exponents of the various types of structure functions 
may help in this regard. This will be attempted in the next section. 

III. ANALYSIS USING EXTENDED SELF- SIMILARITY 

In this section, we shall examine ESS of scalar statistics in the light of the new objects 
defined in terms of their parity. We first apply ESS as is normally done for velocity statistics: 
plot structure functions of all orders against the absolute third-order structure function. We 
shall follow this practice for illustrative purposes, even though the third-order does not 
have a comparably significant meaning for temperature. In Fig. |3(aj| the ESS plot of Si (r) 
versus S\ 3 \ shows a single relative scaling exponent outside the dissipative range. The inset 
shows the local logarithmic slope of the extended self-similarity plot which is calculated via 
D2,\3\ = dlog(S , 2 (r))/dlog(S'|3|(r)). However, this same signal has two scaling regions with 
distinct exponents when plotted against r (corresponding to the passive range at small scales 
and the convective range at large scales, see top panel of Fig. 1). This two-exponent scaling 
behavior is masked by ESS. This is essentially so even if we plot the sign-antisymmetric 
structure function of the second-order against S3 (which is also sign-antisymmetric, see Sec. 
2), as seen in Fig. |3(b)| 

Let us now compare objects of the same order, but different sign-symmetry. In Fig. |4(a)[ 
we show a plot of Sz(r) which is odd-parity (made of the antisymmetric combination of 
the PDF functions) versus SWr), which is the corresponding even-parity object of the 
same order. The two objects scale identically in the smallest (dissipative/diffusive) scales 
with relative exponent of unity, transitioning into a region of relative exponent 1.27. This 
corresponds to the inertial range scaling, as we already know from Fig. |21 Past this range, 
the relative scaling exponent drops back to ~ 1. In Fig. |4(b)"[ we repeat the ESS comparison 
for the second-order. Instead of plotting the second-order object against another order, 
we choose to compare the even- and odd-parity manifestations S2, s (r) = S2 (r) and S2 >a (r), 
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respectively. Once again, it is seen that there are two regions of scaling beyond the dissipative 
range. The small-scale range has a relative scaling exponent of about 1 while the large-scale 
range has a slope of 1.42. These exponents are as expected from the direct scaling analysis 
of fl. 

Finally, when we compare objects of different order and different parity as shown in Fig.|SJ 
the dual scaling range feature persists. Figure [5(aJ| shows S-2 )S = ^(r) versus S^ a = S^{r) 
while Fig. |5(b)| shows S^a vs. S^ s = S^^r). Again, we recover two separate scaling ranges, 
one inertial and one convective. 

From these examples, we can now make the following general statement. In a process with 
two distinct scaling ranges, ESS cannot distinguish between them whenever a comparison 
is made between two like-parity objects. However, if ESS compares objects of any order, 
but of opposite parity, a second scaling (if it exists) can be recovered. This is our major 
qualitative conclusion. This kind of comparison has been made possible by the introduction 
of the odd-symmetry, even-order objects S n>a (r) (n even) — which, to our knowledge, have 
not been considered before. 



IV. PARITY OF THE NEW MOMENTS FROM THE VIEWPOINT OF SO (3) 
DECOMPOSITION 



We shall now discuss the above property of ESS from another perspective that contains 
the correct description of the parity property. We consider an analytical expression for 
structure functions with two distinct scaling ranges, separated by a crossover scale l c . Since 
we do not know the proper analytical expression for the structure functions, we use an 



ll| in this direction has been extended 



empirical interpolation formula. Batchelor's attempt 
variously, in particular in Refs. for structure functions of all orders. Here, we shall 

extend the specific form proposed in Ref. Q| to the nth order as 



SJr) 



9[ h 



-| a n -fi n 



with the function 



9\ x ) 



(18) 



(19) 



(1 + x K )V K ' 

where the assumed scaling exponents a n hold for the passive range, r <C l c , and j3 n for the 
convective range, r 3> Z c , and a n > j3 n ; g(x) is a dimensionless function that is monotonically 



decreasing, and the exponent k > 1 determines the width of the crossover from one scaling 
range to the other. 

We make use of the spatial parity of the signed structure functions by noting that, in 
homogeneous turbulence, S n>s (r) = S\ n \(r) is even-parity with respect to reflection under 
r — > — r for all n; the odd-order normal structure functions S2k+i, a ( r ) = 5 , 2fe+i( r ) are odd- 
parity; the newly defined S 2 k, a , even-order sign-antisymmetric statistics are odd-parity as 
well. So, in the homogeneous case the sign-symmetry with respect to the increment is the 
same as the sign-symmetry with respect to the spatial orientation of r (i.e. parity) and this 
allows for the following decomposition of the objects defined in Sec. II. 

The recently developed SO (3) group decomposition jl^ applies conveniently to the sym- 
metric and antisymmetric structure functions. For the scalar case, the basis functions are 
spherical harmonics Yi >m (r) Q]. The even or odd spatial parity is carried by the angular 
dependence in the spherical harmonics, in effect by the index I, as 

= (-!)%*(*) (20) 

Perfectly isotropic objects would contain an I = sector only, i.e., no angular dependence 
would be present. Sign-symmetric functions are composed of only even sectors while sign- 
antisymmetric functions are composed only of odd sectors. Thus, 

fi n , s (r, 9, 0) = Sg(r) + Sg)(r, 9, <j>) + . . . , (21) 

°° 2k (2k) 

k=0 m=-2k 

where the superscripts (0), (2), ... denote the even-parity contributions allowed, and 

S n>a (r, 6, 0) = fig) (r, 9, 0) + (r, 9, </>) + ..., (22) 

00 2k+1 (2k+l) 

Yl A 2 k+i, m r Cn Y 2k+1>m (9,(j)), 

k=0 m=-2k-l 

where the superscripts (1), (3), ... denote the odd-parity contributions allowed. We can 
substitute the algebraic scaling form of (fTHj) for each scaling term in and to obtain 

(2*) 

S n>a {r, 9, 4>) = ^ A 2k,m 



k,m 

9 



x 



(2k)_ R (2k) 



r 

x Y 2kjm {9,(f)), (23) 
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and 

S n>a {r, Q,4>) = Y1 Mk+u 



(2fc+l) 



fc,m 

(2fc+l) o(2fe+l) 

X 



r 



x F 2fc+1>m (0,0). (24) 

The finding that ESS yields the same relative exponents, when comparisons are made of 
even-parity statistics, means that for all sectors and for any pair of orders (p, q) we have 



(2k) R (2k) 
VLq fjq 



(25) 



Equivalently, for p < q, it follows that 



a (2fe) fj(2k) 

*? k) nfi (i + *P } ) = /f fc) n-/(i + *? fc) ) ' (26) 

where the quantities 

(*) _ (2k) 

f= w V? ■ < 27 > 

are the relative increments between successive orders. The constraint (|25jl means that the 
relative increments in both convex functions, a p (p) as well as (3 p (p), are equal. This was 
already shown in Ref. where physically different turbulent states such as the Kolmogorov 
turbulence in three dimensions, thermal convection as well as magnetohydrodynamic turbu- 
lence showed the same relative ESS exponents. We suppose, therefore, that as long as the 
different regimes belong to the same "universality class" — i.e., if they have the same rela- 
tive increments 5^ for all orders n — ESS will simply mask the presence of distinct scaling 
regimes. 

The argument for the odd-parity ESS comparisons is similar, now with the constraint 

(2fc+l) o(2*+l) 

(2fc+l) n(2k+l) ■ 

Uiq fJq 

One remains within a universality class when plotting S Pta (r) as functions of S q>a (r), but it 
can be expected that the relative exponents will be different from those for the even parity 
case. 
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We now turn to the case where S PiS (r) are compared to S qja (r), which are objects with 
different symmetries. When proceeding as in (|25|h ()26j) . and (|2Tj) we obtain 

4 2k) 



0,(3*0 
P 



4 2fc+1) nS (l + ef k+l) ) $ k+l) in= p (l + e? k+1) ) ' 
For insensitivity to different scaling regimes, the following must hold true: 



(29) 



a< 2 *> $ 2fe > 



(2fc+l) o(2fc+l) 



(30) 



Recall that the superscript (2k) corresponds to even parity objects while (2k + 1) corre- 
sponds to odd parity. We expect that several mechanisms, such as the breaking of reflection 
symmetry by a mean gradient for scales r < l c or by buoyancy effects for scales r > l c , will 
lead to different ratios on the left and right hand side of the Eq. (|30|). respectively, thus 
leading to a "discontinuity" in scaling. While this discussion cannot be a proof of our main 
observation, it shows that the two scaling regimes belong to different universality classes 
when different parities are involved. Further, it provides hints for what kind of quantities 
have to be investigated in order to shed more light on the issue, in experiments as well as 
in numerical simulations. 



V. CONCLUSIONS 



We have analyzed moments of temperature increments with respect to their sign- 
symmetry properties, and defined sign-symmetric and sign-antisymmetric components for 
both even and odd-order structure functions. ESS analysis of all combinations of symmetry 
and order of structure functions indicates that this technique masks the convective scaling 
regime. Only if objects of opposite parity are compared can one recover the distinct scalings. 

We have presented a model for how such scaling regimes behave in the ESS analysis using 
even and odd orders of a spherical harmonic expansion. Unfortunately, this model cannot be 
taken to its logical conclusion because the various numerical coefficients cannot be obtained 
with any certainty from the existing data. Furthermore, a dependence on the strength of 
the mean temperature gradient might modify some of our results as well — e.g., differences 
in the exponents of even order moments in comparison to odd order moments. 

One purpose here has been to point out the possible pitfalls in using ESS without proper 
consideration of the symmetry properties of the statistics being compared. This observation 
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does not detract from the merits of the method. Indeed, ESS has proved to be a very useful 
tool in extracting exponents when the scaling range with respect to r is short. 
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FIG. 1: Normal structure function of the temperature fluctuations. Top figure: S^r) is shown for 
height h = 1.75 m above the ground. Shown by the vertical arrow on the abscissa is the integral 
scale Lq for the temperature fluctuations. Bottom figure: S^r) as a function of S^r). The data 
corresponding to the top figure are marked as circles (multiplied, for convenience, by a factor of 
10). A and B, which have the same slope in the bottom figure, correspond to the two scaling 
ranges with slopes of 0.65 and 0.44, respectively, in the top figure. The crosses in the bottom figure 
correspond to another similar situation where A and B, again, have two distinctly different slopes 
for normal structure functions. Both figures are reproduced from Ref. [J. 
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(a) Logarithmic local slopes of the various (b) Logarithmic local slopes of the various 

structure functions for n = 2 structure functions for n = 3. 

FIG. 2: Local slopes in the log-log graphs of four types of structure functions: positive (P n ) 
and negative (N n ) parts, the sum of positive and negative parts (same as the absolute moment), 
S\ n \ = Pn + N n = S njS , and the negative minus the positive part, \S n \ = N n — P n = \S nja \- 
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(a) Loglog plot of S 2 vs. S\ 3 \ (b) Loglog plot of S 2 , a vs. S 3 . a 

FIG. 3: ESS comparisons of structure functions of different order and same sign-symmetry. Insets 
show the (logarithmic) local slope in these ESS coordinates. In this comparison, only one scaling 
exponent is seen past the dissipation range. The conditions for this and other figures to follow are 
the same as for Fig. 1. 
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(a) Loglog plot of S3 vs. 5| 3 | 



(b) Loglog plot of S 2 , a vs. S 2 ,, 



FIG. 4: ESS comparisons of structure functions of the same order and different sign-symmetry. 
Inset shows the (logarithmic) local slope in these ESS coordinates. For this comparison, the two 
different scaling regimes are retreived. Vertical bars are inserted in order to show regions of constant 
local slope. 
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S 3,a S 3,s 



(a) Loglog plot of S 2 , s vs. 5 3ia (b) Loglog plot of S 2 , a vs. 5 3jS 

FIG. 5: ESS comparisons of structure functions of different order and different sign-symmetry. 
The insets show the (logarithmic) local slope in these ESS coordinates. Once again, two different 
scaling regimes are retrieved. Vertical bars are inserted in order to show regions of constant local 
slope. 
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